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We describe a variety of intriguing mode-coupling effects which can occur in a 
confined Bose-Einstein condensed system at finite temperature. These arise from 
strong interactions between a condensate fluctuation and resonances of the thermal 
cloud yielding strongly non-linear behaviour. We show how these processes can 
be affected by altering the aspect ratio of the trap, thereby changing the relevant 
mode-matching conditions. We illustrate how direct driving of the thermal cloud 
can lead to significant shifts in the excitation spectrum for a number of modes and 
provide further experimental scenarios in which the dramatic behaviour observed 
for the m = mode at J1LA (Jin et al. 1997) can be repeated. Our theoretical 
description is based on a successful second-order finite-temperature quantum field 
theory which includes the full coupled dynamics of the condensate and thermal 
cloud and all relevant finite-size effects. 
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1. Introduction 

Bose-Einstein condensation (BEC) in dilute atomic gases offers a virtually unique 
laboratory for the investigation and experimental testing of thermal quantum field 
theories (QFT). The ability to perform both a priori calculations and precision 
experiments for direct comparison has led to a deepening of our understanding of 
the important processes which must be included in a successful finite-temperature 
QFT. The pioneering measurements of condensate excitation frequencies and de- 
cay rates at the Joint Institute for Laboratory Astrophysics (JILA) have provided 
the most stringent tests of such theories to-date (Jin et al. 1997) and efforts to 
explain these results have led to the development of a second-order QFT which 
simultaneously includes both the dynamic coupling of the condensed and uncon- 
densed atoms and the effects of strong interactions and finite-size (Morgan et al. 
2003, Morgan 2004). Of particular importance is the fact that the theory explicitly 
includes the effect of direct driving of the thermal cloud by the perturbation which 
generates the excitations. This allows the condensate to be excited indirectly via 
the non-condensed atoms as well as directly via the perturbation, a process which 
can be very important at finite temperature when there is a large non-condensate 
fraction. The final agreement between the theoretical predictions and the JILA 
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measurements is excellent, with no adjustable parameters required in the theory 
(Morgan et al. 2003). This was the first QFT calculation for confined BECs which 
systematically included the phonon character of the low lying states, all relevant 
Landau and Bcliaev processes and the effects of finite system size. 

In this paper we propose new spectroscopic experiments to further test this 
thermal QFT and to explore the intriguing mode-coupling effects which can occur 
in a trapped condensate at finite temperature. Initially we consider a condensate 
confined in a spherically symmetric harmonic trap. This geometry has been pre- 
viously considered (Rusch et al. 2000, but in the absence of direct driving of the 
thermal cloud and by the Pisa group Liu et al. 2004) using a different formalism 
based upon the random phase approximation. The importance of this geometry 
is in the mode matching characteristics of the excitation spectrum. In a spherical 
trap, a condensate has a breathing mode excitation with zero angular momentum 
(I = 0) and a frequency which lies just above a strong resonance of the thermal 
cloud. As a consequence, there is a strong coupling between the two with the result 
that the mode energy is almost impervious to temperature, increasing slightly as 
the system heats up in contrast with the substantial downward shifts predicted 
for all other modes studied to date and observed in experiment (Jin et al. 1997, 
Morgan et al. 2003, Hutchinson et al. 1997, Dodd et al. 1998, Giorgini 1998, Mor- 
ago et al. 2001). In addition, at temperatures approaching the critical temperature, 
where the thermal cloud dynamics can become dominant, indirect excitation of the 
condensate via the thermal cloud as intermediary can lead to a dramatic shift in 
the I = excitation frequency. This is in direct analogy to the similar behaviour of 
the lowest-energy m = mode observed in the JILA TOP trap experiment which 
provoked much theoretical discussion before its final resolution (Morgan et al. 2003, 
Jackson & Zaremba 2002). 

Since this strong coupling is a consequence of mode matching between conden- 
sate and thermal modes it should be possible to adjust this behaviour by tuning the 
condensate excitation frequency closer to that of the thermal cloud. This can be 
achieved in practice by altering the anisotropy of the trapping potential (Hodby et 
al. 2001). In particular, guided by the zero-temperature mode spectrum calculated 
by Hutchinson and Zaremba (1998), we identify a resonance between a low-lying 
m = condensate mode and the thermal cloud with a frequency u> = 2oj r for a trap 
with an aspect ratio of u> z /u> r = 0.764, where u> r and u> z are radial and axial trap 
frequencies respectively. We apply our formalism in this geometry and demonstrate 
the interesting phenomena which can be expected in this case. 

Observation of the predicted effects in a spherical trap, and particularly in 
this resonant case, would provide additional confirmation of the validity of our 
theoretical description and intriguing further examples of the coupled dynamics of 
condensed and uncondensed atoms which can occur in trapped Bose gases at finite 
temperature. 
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2. Theory 

Briefly, (Morgan et al. 2003, Morgan 2004) we begin from the generalized Gross- 
Pitaevskii equation (GPE) for a condensate of Nq atoms 

ih^- = [F SiD + p(r,i)-A(;) + iVo(;)[/o|$| 2 ]$ 

+ 2U n{r, i)$ + U rh(r, t)$* - /(r, t), (2.1) 

where H sp = — — h V^^r) is the single particle Hamiltonian, containing the 
kinetic energy and the static trap potential, V ex t(r). We take V ex t( T ) to be the 
usual cylindrically symmetric harmonic potential V r ex t(r) = ^m(uj 2 r 2 + uj 2 z 2 ) with 
r 2 = x 2 + y 2 . The other terms in Eq. (2.1) are as follows: P(r, t) is a time-dependent 
external perturbation which generates excitations, X(t) is a scalar which controls 
the condensate phase, n(r, t) is the non-condensate density, m(r, t) is the so-called 
anomalous average and /(r, t) is a term resulting from particle number conservation. 
Interactions are described by the usual contact potential with interaction strength 
Uq = Ai:h 2 a s /m where a s is the s-wave scattering length. This means that our 
definition of m(r,i) is subject to the usual ultra-violet renormalization. 

In the static case, the GPE has a time-independent solution which satisfies 

[H sp -X + N U \$(r)\ 2 + 2U ri(r)Mr) 

+ E/bm(r)S*(r)-/(r)=0, (2.2) 

where A is the condensate eigenvalue and is approximately equal to the chemical 
potential (Morgan 2000). The associated quasiparticle excitations obey a set of 
coupled Bogoliubov-de Gennes equations which we solve for the static quasiparticle 
wavefunctions {«i(r), Vi(r)} and energies {e^}. We treat the thermal cloud as a 
non-interacting gas of quasiparticlcs coupled to the condensate mean-field with the 
result that all thermal cloud quantities can be calculated directly in terms of these 
quasiparticle solutions. 

We now consider the effects of the external perturbation P(r,i), allowing all 
quantities to develop small time-dependent oscillations about their static values; 
<&(r,i) = <f>(r) + <5$(r,i), h(r,t) = n(r) + 5n(r,t) etc. Substituting into the time- 
dependent generalized GPE, Eq (2.1), and linearising we obtain the equation of 
motion for the condensate fluctuation <5$(r, t). This is solved by combining it with 
its complex conjugate, Fourier transforming and expanding in the static quasipar- 
ticle basis 

The bi(u>) expansion coefficients are directly related to the experimentally measured 
condensate density fluctuations, given by Sn c = 5(N \$>\ 2 ). 

Assuming that only a single condensate mode "p" is excited, the expansion 
coefficients are of the form 

b p (u) = P p() (uj)G p (uj + ij), (2.4) 

where P p o(lo) is the matrix element for the generation of the excitation from the 
condensate and ij is a small imaginary part in the frequency (see later) . The resol- 
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vent G p (tu) is defined as 

= * 1 v ( V ( 2 ' 5 ) 

nu — e p — ±j p (uj) 

where H p (u>) is a frequency-dependent self-energy given by 

E p (u J )=AE p f+AE^(co). (2-6) 

E p (w) contains both static (S) and dynamic (D) contributions to the energy. The 
static part, AE PP , describes interactions between the condensate and the equi- 
librium non-condensate mean fields, while the frequency-dependent dynamic part 
AE P ^ (tu) describes the driving of the thermal cloud by the condensate fluctuation 
and the subsequent back-action. 

So far we have assumed that the external perturbation only directly excites the 
condensate. All the dynamics of the thermal gas is therefore due to the subsequent 
action of the condensate on the uncondensed cloud. A second possibility exists, 
however, in which the external perturbation also directly drives the non-condensate. 
Including this process, we find that the expansion coefficients are given by 

bp(w) = P p0 {lu)R p {uj + ij), (2.7) 

where the new resolvent R p (lu) is defined by 



PpQ M 



G p (u). (2.8) 



Again there are both static and dynamic contributions present in this expression. 
The static term AP^p is due to the effect of changes in the static condensate shape 
on the excitation matrix element, while the more interesting dynamic term is due to 
direct excitation of thermal cloud fluctuations by the external potential and their 
subsequent coupling to the condensate. 

Both the dynamic terms AE pp \lu) and AP p q(u>) are calculated via a sum 
over all the Landau and Beliaev processes which can occur in the thermal cloud. 
These are resonant whenever certain energy matching criteria are satisfied and the 
parameter 7 is therefore required to keep the expressions finite at the corresponding 
poles. Physically, 7 describes the finite resolution of the experiment and typical 
values are of order a few times 10 -2 of the trapping frequency. In this paper we 
have taken 7 = 0.05aj r but our results are not sensitive to its precise value within 
the experimentally relevant range. 



3. Spherical Geometry 

Figure 1 shows results for the energies of three quasiparticle excitations in a spheri- 
cal trap as a function of reduced temperature. The condensate contains 6000 atoms 
and the trap frequency is uj r — 2ir x 100Hz. Results are shown from our second 
order QFT including and excluding the effect of direct thermal driving, as well 
as from a perturbative HFB-Popov calculation in which the thermal cloud is held 
static (Griffin 1996, Hutchinson et al. 1997, Hutchinson et al. 2000). The external 
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Figure 1. Excitation energies as a function of reduced temperature t = T/T° for the 
lowest-energy (n = 1, I = ), (n = 0, I = 3) and (n = 0, / = 2) quasiparticle modes, 
from highest energy down. The diamonds and circles exclude (G) and include (R) thermal 
driving, respectively. Squares give the mean-value of the self-energy at the poles of the 
resolvent, the dashed lines are the unperturbed Bogoliubov energies and the solid lines 
are the results of a perturbative implementation of the HFB-Popov theory. Energies are 
in units of the radial trapping frequency. 

perturbation was assumed to be a low-order polynomial multiplied by the spherical 
harmonic with the appropriate symmetry for each mode. For the I = and I = 2 
modes it therefore had a quadratic dependence on the co-ordinates, while for the 
I = 3 mode it was a cubic function. 

A number of features of Fig. 1 arc of interest and demonstrate the importance of 
the dynamic coupling between the condensate and non-condensate. In general, the 
second order theory predicts markedly greater shifts than the HFB-Popov theory. 
For the I = 3 and I = 2 modes, this shift is negative at high temperature and direct 
driving of the non-condensate has little effect. In contrast, the I = mode is almost 
unaffected by temperature, shifting slightly upwards as the temperature rises, in 
contrast to the HFB-Popov prediction which shows a substantial downwards shift. 
For this mode, there is also a dramatic difference between the results obtained 
with and without direct driving of the thermal cloud. If this process is included, 
the response shifts towards the non-interacting gas value of 2u) r for temperatures 
greater than about 0.7T C °, where Tj? is the ideal gas critical temperature. This 
is reminiscent of the similar shift seen for the m = mode in the TOP trap of 
the JILA experiment and the physical explanation is the same in this case; the 
condensate excitation is usually excited directly by the external perturbation but 
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at high temperatures it can also be excited indirectly via the thermal cloud. If the 
thermal cloud response is particularly large at the frequencies of interest and if the 
coupling to the condensate is sufficiently strong, this second excitation method can 
become dominant and the condensate response becomes peaked around that of the 
non-condensate. 

We can gain some insight into the effects of the coupling between the condensed 
and uncondensed atoms by considering a simple model in which the thermal cloud 
response is characterized by a single resonance at a frequency ujq. If the condensate 
has a Bogoliubov energy of E p = hu> p and the (temperature-dependent) coupling 
strength to the thermal cloud is g 2 then the resolvent G p (<jj) has the form 

ti8 p (u>) = ; -5—. (3-1) 

to + ij — oj v 

In the limit g — 0, this is the usual lorentzian with a single peak at the Bogoliubov 
frequency, to = uj p . For finite g and in the limit 7 — > 0, Q p has two poles at u> = 
(u} p + ujo)/2±y / (ui p — wo) 2 /4 + g 2 . Thus as g increases the poles move further apart 
and in particular the original pole at u> ~ uj p is shifted away from the resonance at 
wo- This is the usual level repulsion familiar from second order perturbation theory. 

If we include the effect of direct driving of the non-condensate with coupling 
strength gp, the resolvent Q p (to) is replaced by JZ p (u)), which has the form 

I _|_ spg 

hHp(w) = " +t7 '"° 2 — . (3.2) 

U) + i7 — UJn 7-9- 

If gp is small and/or the frequency difference lo p — ujq is large, this does not differ 
greatly from Q p (u>) for frequencies around lo p . However, if gp is large and u> p — luq 
is small, the additional factor in the numerator can dominate the response which 
becomes again a simple lorentzian but now centred on the thermal cloud resonance 
u)q. This cross-over behaviour can clearly be seen in Figs. 2a and 2b where we 
plot the response functions Q p {to) and 1Z p (lu) as a function of frequency for various 
reduced temperatures. The strong response at uj = 2to r is visible in 1Z p (oj) at 0.9T® 
but is absent in Q p (lu) which actually has a local minimum in this region. At this 
temperature and frequency, the excitation of the condensate via the thermal cloud 
is more than four times as large as the direct effect of the perturbation. 

The frequency-dependence of the real and imaginary parts of the self-energy 
E p (w) for this mode are shown in Fig. 2c) at a temperature T = 0.9T°. Although 
this clearly has other structure, its form is rather close to that assumed in our 
analysis above, consisting of a (negative) constant and a single large resonance at 
a frequency u) = 2lo,.. The value of the negative constant for this case is of order 
— 0A3huj r and without the resonance structure the mode would shift downwards 
by about this amount. Comparison with Fig. 1 shows that such a shift is of the 
same order as the downwards shifts observed for the / = 3 and / = 2 modes 
at this temperature. The reason the I = mode is instead almost impervious to 
temperature in the absence of direct thermal driving is therefore the result of strong 
level repulsion from the thermal cloud resonance at lo — 2u> r . The upward shift this 
produces almost exactly counterbalances the downward shift we would otherwise 
obtain from the cumulative effect of all other modes in the problem. That this 
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Figure 2. Response functions and self-energies for the (n = 1,1 = 0) mode in a spherical 
trap, a) Q p (u>) at t = 0, 0.75, 0.9 (solid, dotted, dashed (x2) respectively), b) 1Z p (lj) at t 
= 0, 0.75, 0.9 (solid, dotted, dashed (x2) respectively), c) Real (dashed) and imaginary 
(solid) parts of the self-energy £ p (cj) at t = 0.9. 



interpretation is correct can also be seen from the appearance of a second peak in 
the response function Q p (to) at to ~ 1.7 at t = 0.9 shown in Fig. 2a. Comparing 
with our simple model above, we see that this comes from the second pole which 
originates at the thermal cloud resonance to = 2to r in the absence of interactions 
but is pushed down by level repulsion with the primary condensate resonance at 
lu ~ 2.2hto r . 

The thermal cloud resonance at to = 2to r has its origin in a large number of 
Landau processes which correspond to removing an atom from a single-particle 
level with quantum numbers (m, h) and placing it in a level with quantum numbers 
(n2,h)- Since the energies of an ideal gas in a spherical trap are given by e n i = 
hto r (2n + l + 3/2), these processes are resonant at frequencies to 12 — to r (2(n 2 — ni) + 
h — h), so a t high temperature the non-condensate has a large response at integer 
multiples of the trap frequency. In particular, there is a large response around to = 
2to r corresponding to processes where 1% — h is even. Since the perturbation required 
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to excite a Bogoliubov mode with angular momentum l p involves the spherical 
harmonic Y l p (8, </>), we have the parity selection rule l p + h + h = even. Thus, if 
the condensate has a quasiparticle excitation with even angular momentum close in 
energy to an even integer multiple of the trap frequency we can expect its dynamics 
to be strongly influenced by the thermal cloud. This situation indeed occurs in a 
spherical trap where the lowest I = mode has an energy of order 2.2h\o r . The 
I = 3 mode shows no such influence of the thermal cloud, despite having an energy 
significantly closer to 2hw r , because it has odd parity and is therefore decoupled 
from the thermal cloud resonances in this region. 

In a TOP trap a similar situation arises because there is a mode with even 
z-parity and zero angular momentum about the z-axis (m = 0) with a Bogoliubov 
energy of order 1.8Sw r . This mode is therefore strongly coupled to Landau reso- 
nances with frequencies around 2lo t and the same effect of direct driving of the 
non-condensate is observed. In this case, however, the condensate mode lies below 
the thermal cloud resonance with the result that the effect of thermal cloud driving 
is to shift its frequency up towards the non-interacting value, in contrast to the 
downwards shift this process produces for the I = mode in a spherical trap. This 
difference also shows up in the relative phase of the condensate-non-condensate 
oscillations. When the condensate drives the thermal cloud these oscillations are 
generally out-of-phase at high temperature, but in the reverse situation the relative 
phase depends on whether the thermal cloud resonance is above or below the unper- 
turbed condensate mode. For the m = mode at JILA this resonance is above the 
primary condensate energy and the growing importance of direct thermal driving is 
associated with a transition to in-phasc oscillations of the two clouds. In contrast, 
in the spherical trap the thermal cloud resonance is below the original condensate 
mode and we find that the oscillations remain predominately out-of-phase at the 
highest temperatures. 

We note that whereas our results for the self-energy and resolvent Q p {lo) are 
intrinsic to the system, depending only on the geometry and particle interactions, 
the effects of direct thermal driving are sensitive to the particular form of the 
perturbation used to excite the system and so will vary from one experiment to 
another. This provides an important handle on the effects we have discussed since 
it should be possible to vary the strength of direct thermal driving by modifying the 
form of the perturbation. For example, it should be possible to design an experiment 
to probe our predictions based upon either a localised perturbation of the magnetic 
trap or a Raman scheme. The localised perturbation could be designed to couple 
predominantly to the condensate, which occupies the centre of the trap, rather 
than to the non-condensate which is peaked nearer the edges of the atomic cloud. 
A Raman based scheme could be designed to be resonant with the condensate and 
not the non-condensate, thus directly exciting only the condensate atoms. Such 
experiments would allow the effects of direct thermal driving to be separated from 
other self-energy interactions and give access to the intrinsic response function 

4. Resonant Geometry 

To explore mode matching effects further, we can tune the condensate mode into 
resonance with a thermal cloud excitation frequency by adjusting the aspect ratio 
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of the trap. This is experimentally feasible and has been used to explore nonlinear 
effects such as parametric down-conversion in the past (Hodby et al. 2001). By 
selecting an aspect ratio of ui z /uj r = 0.764 we are able to tune the zero temperature, 
Bogoliubov frequency for one of the low lying m = modes to be precisely 2u> r 
(this mode corresponds to the quantum numbers (rra,p, n) — (0, 0, 2) where m is the 
angular momentum about the z-axis, p is the z-parity (0 is even) and n indicates 
the energy in this subspace). With this geometry we repeat the calculation and 
again see shifts in the excitation frequencies of low-lying modes with increasing 
temperature as displayed in Fig. 3. 



T 




0.1 0.3 0.5 0.7 0.9 

t = T/T° 

c 

Figure 3. Excitation energies for the (m,p,n) — (0,0,2) (lower) and (0,0,3) (upper) 
modes in an anisotropic trap as a function of reduced temperature. Symbols are as in 

Fig. 1- 

It is interesting to note that when we tune the (0, 0, 2) mode into resonance, 
due to a fortuitous anti-crossing behaviour, the (0,0,3) mode is relatively close in 
frequency. We therefore see an additional strong coupling effect for this mode. In 
fact we find that because of the change in matrix elements associated with the 
change in trap geometry, this mode couples more strongly to the thermal cloud 
resonance at uj = 2ui r than the mode tuned to resonance. Its behaviour is therefore 
very similar to that previously observed for the I = mode in the spherical trap 
and plots of the response functions Q p {uj) and 1Z p (uS) and self-energy Yi p (ui) for this 
mode are qualitatively similar to those given in Fig. 2. 

In contrast, the coupling in the self-energy for the (0,0,2) mode is weaker (al- 
though still quite strong) with the result that this mode has the usual strong down- 
wards shift in energy in the absence of direct thermal driving. The mode is nonethc- 
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less strongly affected by direct excitation of the thermal cloud, however, so when 
this process is included we get a strong component to the response around to = 2tu r . 
In contrast with the other cases we have examined, this now causes the response 
to split into two distinct peaks. This is shown in the upper panel of Fig. 4 where 
the resolvents, Q p {lo) and lZ p (u)) are plotted at a temperature of 0.7T®. The two 
peak structure comes from the product of the underlying lorentzian function for 
G p (u) centred near u) = l.9uj r and the thermal cloud resonance at 2.0w r which 
modifies Q p (cu) to TZ v (uo). The product of these two misaligned peaks yields the 
bimodal structure seen in Fig. 4. Note that since this bimodal distribution renders 
a lorentzian fit meaningless, a single frequency cannot be ascribed so this is omitted 
in Fig. 3. 



0.5 




Figure 4. Upper: 7 x \Q p (lu)\ (dashed) and 7 x \1Z p (ijj)\ (solid) vs frequency at reduced 
temperature t — 0.7 for the (0,0,2) mode in a resonant trap. The vertical dashed line shows 
the position of the original unperturbed Bogoliubov mode. Lower: Plot of the condensate 
density response (arbitrary units) in the time domain at t = 0.7 including the effect of 
direct driving of the non-condensate (corresponds to the solid line above). Collapse and 
revival of the condensate oscillations are clearly visible as a result of the split of the 
spectrum into two distinct parts. 

A double-peak response of this form is a signature of collapse and revival between 
coupled modes. This is most clearly seen in the time domain shown in the lower 
panel of Fig. 4. Here we have plotted the condensate response for the (0, 0, 2) mode 
as a function of time for the same temperature as in the upper panel. Initially the 
oscillations decay, but after approximately ten trap periods the amplitude increases 
again as the mode intensity is transfered back from the thermal cloud. We should 
note that the coupling here is to a resonance formed from the cooperative effect 
of many Landau processes. This means that damping could become a significant 



Article submitted to Royal Society 



Novel Dynamical Resonances in Finite- Temperature Bose-Einstein Condensates 11 

issue, although our estimates of the damping rate indicate that this collapse and 
revival should be observable on the time scales of typical experiments. 

The fact that this geometry produces two modes with the same symmetry close 
in energy and both strongly coupled to a thermal cloud resonance means that in 
principle it may be possible to see interesting new effects in this case where energy is 
transfered between the two modes via the thermal cloud as mediator. Such processes 
are included in the QFT we have developed (Morgan 2004) but require a significant 
extension of our numerical methods and are therefore not contained in the present 
calculation. Such effects along with those described in this paper are part of the 
rich dynamical structure which exists in condensed systems at finite temperature 
and which is now accessible to detailed experimental and theoretical study. 

5. Conclusions 

In conclusion, we have presented concrete predictions for the excitation energies of 
low-lying modes of a Bosc-condcnsatc in a spherical trap at finite temperature using 
a successful second order QFT which includes the dynamic coupling of condensed 
and uncondensed atoms and finite size effects. For the lowest breathing mode with 
zero angular momentum, we have observed the effects of strong coupling between 
the condensate and thermal cloud. This leads to a mode energy which is almost 
independent of temperature, shifting only slightly upwards, compared with the 
substantial downward shifts predicted and observed in previous calculations and 
experiments. In addition, at high temperature we observe a strong crossover to the 
non-interacting limit which is a consequence of indirect excitation of the condensate 
via the thermal cloud as intermediary, in direct analogy to the similar phenomenon 
observed in the JILA TOP trap. The observation of these effects would provide fur- 
ther confirmation of the validity of our theoretical approach and of the importance 
of the dynamic coupling between the condensate and thermal cloud in trapped gas 
BECs at finite temperature. 

We have also identified a new regime where, by altering the aspect ratio of the 
trap, it is possible to bring condensate modes into resonance with modes of the 
thermal cloud. Strong mode coupling, especially when the thermal cloud is directly 
driven by the external perturbation, is then observed allowing the observation of 
non-linear phenomena such as collapse and revival between condensate and thermal 
cloud modes. We also note that, due to a fortuitous anti-crossing in the condensate 
mode spectrum, another mode with the same symmetry is also coupled to the ther- 
mal cloud resonance in the vicinity of lu = 2uj r . This should allow for a coupling 
between condensate modes mediated via the non-condensate dynamics which would 
prove a very interesting subject for future theoretical and experimental investiga- 
tion. 

We thank the Royal Society (London), the Marsden Fund of the Royal Society of New 
Zealand and the University of Otago for financial support for this project. 
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